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Monte Carlo calculations have been carried out to determine the absorbed 
dose and dose equivalent for 592-MeV protons incident on a cylindrical phantom 
and for neutrons from 580-MeV proton-Be collisions incident on a semi-infinite 
Phantom. For both configurations, the calculated depth dependence of the absorbed 
dose is in good agreement with experimental data. 


INTRODUCTION 

Nucleon-meson cascade calculations were re- 
cently carried out for monoenergetic neutrons (60 
to 3000 MeV) and protons (U00 to 3000 MeV) normal- 
ly incident on a semi -infinite slab of tissue (ref. 
1). For several incident energies, the calculated 
depth dependence of the absorbed dose was compared 
with available experimental results (ref. 2), 
which were obtained using somewhat different source- 
geometry configurations, and substantial differ- 
ences were found. It was not clear from these com- 
parisons whether the theory was in error or whether 
the differences could be ascribed to the different 
configurations used in the calculations and experi- 
ments . 

The calculations have been carried out for two 
configurations: 592-MeV protons incident on a 

cylindrical phantom, and neutrons from 580-MeV 
proton-Be collisions incident on a semi-infinite 
slab phantom. The calculated depth distributions 
for the absorbed dose are compared with the mea- 
sured distributions of Baarli and Goebel (ref. 2) 
for these same configurations. The depth distribu- 
tions for the dose equivalent have also been calcu- 
lated . 

The method of calculation used is summarized 
in the next section, and the results are presented 
and discussed in the last section. 

CALCULATIONAL METHOD 

The method of calculation is the same as that 
described previously (ref. l), so only a brief ac- 
count of the c alculat i onal method will be given 
here. The calculations were carried out using the 


Monte Carlo code NMTC (ref. 3). This code takes in- 
to account charged-particle energy loss due to ioni- 
zation and excitation of atomic electrons, elastic 
q.nd nonelastic nucleon-nucleus and pion-nucleus col- 
lisions, pion and muon decay in flight and at rest, 
and negative-pion capture at rest. The electron- 
photon cascade resulting from the decay of neutral 
pions and the electrons and positrons from muon de- 
cay are included in an approximate manner (ref. l). 
At energies above 15 MeV for nucleons and 2.2 MeV 
for charged pions, the energy, direction, and mul- 
tiplicity of particles produced in nucleon-nucleus 
and pion-nucleus collisions are calculated using 
the intranuclear-cascade-evaporation model of nuc- 
lear reactions (refs. U, 5)* Proton- nucleus col- 
lisions below 15 MeV and pion-nucleus collisions be- 
low 2.2 MeV (except for the capture of negative pi- 
ons at rest) are neglected. Neutron collisions be- 
low 15 MeV are treated using experimental cross- 
section data (ref. 3) in conjunction with the evap- 
oration model (ref. 5) to determine particle pro- 
duction from neutron- nucleus nonelastic collisions. 

A detailed description of the method of calculation 
is given elsewhere (ref. l). 

The calculations Rave been carried out using 
tissue for the phantom composition. The following 
concentrations were used for the tissue (in at. %) : 

H, 63.3; 0, 25.8; C, 9o ; N, l.t. A tissue density 
of 1.0 g/cm 3 was used. 

RESULTS AND DISCUSSION 
Incident Protons 

The geometry for the case of 592-MeV incident 
protons is shown in figure 1. The protons are in- 
cident normally on one end of a cylindrical phantom 
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200 cm in depth by 50 cm in diameter. The experi- 
mental spatial variation of the incident proton in- 
tensity (ref. 2) was approximated in the calcula- 
tions as a Gaussian distribution with a full-width- 
at-half-maximum value of 30 cm. 
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FIGURE 1.— Configuration for incident protons. 


The depth dependence of the absorbed dose is 
shown in figure 2. The experimental absorbed dose 
is not given in absolute units, and thus only the 
shapes of the calculated and experimental distri- 
butions can be compared. Also shown for comparison 
in figure 2 is the absorbed dose calculated for a 
cylinder of infinite radius. This case corresponds 
to the configuration considered in the previous 
calculations (ref. l), i.e., an infinitely broad 
beam incident on a semi-infinite slab. As shown in 
figure 2, the shape of the calculated distribution 
is in very good agreement with the experimental 
distribution when the cylindrical geometry used in 
experiment is simulated in the calculations. There- 
fore, for a meaningful comparison with the experi- 
mental absorbed-dose distribution it is necessary 
that the finite radius of the phantom be taken in- 
to account in the calculations. 
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FIGURE 2.— Depth dependence of absorbed dose for 
592-MeV incident protons. 

The contributions to the absorbed dose and 
dose equivalent from various kinds of particles 
have also been calculated, and these are shown in 
figures 3 and k for the finite-radius geometry. 

The total dose is obtained by adding all of the 
contributions shown in the figures. The histogram 
labeled ’’primary ionization” in figure 3 gives the 
dose from the excitation and ionization of atomic 
electrons by those incident protons which have not 
undergone nuclear collision. The histogram labeled 
’’secondary protons" gives the absorbed dose from 
the excitation and ionization of atomic electrons 
by protons produced from nonelastic nucleon-nucleus 
and pion-nucleus collisions and from the elastic 
collisions of nucleons and pions with hydrogen nuc- 
lei. The histogram labeled "heavy nuclei" gives 
the absorbed dose from particles with mass number 
greater than one produced from nonelastic nucleon- 
nucleus and pion-nucleus collisions and the ab- 
sorbed dose from the recoiling nuclei produced 
from elastic neutron-nucleus collisions and from 
nonelastic nucleon-nucleus and pion-nucleus col- 
lisions. The histogram labeled "charged pions" 
gives the absorbed dose from the excitation and 
ionization of atomic electrons by both positively 

and negatively charged pions produced from nucleon- 
nucleus and pion-nucleus nonelastic collisions. 

The histogram labeled "photons from neutral pions" 
gives the absorbed dose from the electron-photon 
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cascade produced by the photons which arise from 
the decay of neutral pions. The histogram labeled 
"electrons, positrons, and photons" gives the ab- 
sorbed dose from the electrons and positrons pro- 
duced by muon decay and the absorbed dose from the 
photons produced by nucleon- nucleus and pion-nucleus 
nonelastic collisions. The histogram labeled 
"muons" gives the absorbed dose from the excitation 
and ionization of atomic electrons by both posi- 
tively and negatively charged mu-mesons. 



FIGURE 3.— Contribution of various particles to the 
absorbed dose for 592-MeV protons incident on a 
cylindrical phantom. 

The histograms in figure k have meanings simi- 
lar to those of figure 3 but correspond to the dose 
equivalent from various kinds of secondary particles. 
The calculation of the dose equivalent was carried 
out taking the quality factor to be a function of 
the linear energy transfer as in the previous cal- 
culations. In the case of protons, the damage 
curve given in reference 6, which is based on the 
recommendations of the National Committee on Radi- 
ation Protection and Measurements, was used. In 
the case of charged pions and muons, the quality 
factor as a function of linear energy transfer was 
taken to be the same as that for protons , and dam- 
age curves for charged pions and muons, constructed 
in a manner similar to the proton damage curve 
given in reference 6, were used. A quality factor 
of 20 was assigned to the energy deposited by all 
heavy nuclei and a quality factor of unity was 


assigned to the energy deposited by electrons, posi- 
trons , and photons . 



FIGURE Contribution of various particles to the 
dose equivalent for 592-MeV protons incident on a 
cylindrical phantom. 

The contributions to the absorbed dose and dose 
equivalent by various kinds of particles have also 
been calculated for the case of a cylinder with in- 
finite radius, and these results are given in 
reference 7. 

Incident Neutrons 

The experimental arrangement used by Baarli 
and Goebel (ref. 2) to measure the absorbed dose 
from incident neutrons is shown schematically in 
figure 5 • The neutrons produced from proton- Be 
interactions at angles of 18 ° and 56° with respect 
to the proton beam were directed toward an absorb- 
er of tissue-like material, and the depth depen- 
dence of the absorbed dose was measured. The 
neutron intensity on the absorber was approximately 
uniform over a circular area 15 cm in radius. The 
energy spectrum of the incident neutrons was not 
measured. A "nominal" energy for the neutron beam 
was calculated by Baarli and Goebel assuming only 
elastic collisions in the target. These nominal 
energies are 525 MeV and l80 MeV for scattering 
angles of 18 ° and 56°, respectively (ref. 2). In 
the previous calculations for the depth dependence 
of the absorbed dose (ref. l), all of the incident 
neutrons were assumed to be monoenergetic at ener- 
gies of 525 or 180 MeV. 
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FIGURE 5. — Configuration for incident neutrons. 

To obtain an estimate of the incident neutron 
spectrum, the energy distribution of neutrons from 
580-MeV protons incident on Be was calculated using 
the intranuclear-cascade model (ref. U). The re- 
sults, in terms of the differential neutron pro- 
duction cross section, are* shown in figure 6. In 
calculating the neutron production spectrum, the 
production averaged over the angular intervals from 
15. 5° to 20.5° and from 51° to 6l° were used to 
represent the production at 18 and 5 6 , respec- 
tively. It is evident from figure 6 that repre- 
senting these spectra by monoenergetic neutrons at 
180 MeV and 525 MeV, as was done in the previous 
calculations (ref. l), is a gross oversimplifica- 
tion . 


The depth dependence of the absorbed dose for 
the case of incident neutrons is shown in figures 
T and 8. The histograms labeled ”l80-MeV spectrum” 
and "525-MeV spectrum” were calculated using the 
incident neutron spectra shown in figure 6. The 
histograms labeled "l80 MeV” and "525 MeV” are 
taken from the previous calculation (ref. l) for 
monoenergetic incident neutrons • Since the ex- 
perimental results are not reported in absolute 
units, they have been normalized to the present 
calculated results at a depth of 8 cm. The mea- 
sured and calculated depth dependences of the ab- 
sorbed dose are in good agreement when the energy 
distribution of the incident neutrons is taken in- 
to account in the calculations. 

The results in figures 7 and 8 have been nor- 
malized on a per-inci dent-neutron basis. The re- 
sults could have been normalized per proton-Be 
interaction. Letting N be the number of neutrons 
incident on the absorber per proton-Be interaction. 


N = 


Aft T 

a inel o 


J°° dE 


A z a 

AEAft 


where (A 2 a/AEAft) is given in figure 6, o inel is 
the Be inelastic cross section for protons at 
580 MeV (= 210.7 mb, from the intranuclear-cascade 
calculations), and Afl* is the solid angle for neu- 
trons from the target hitting the absorber, which 
can be obtained from the geometry of the experiment 
[Aft* = tt(15) 2 /(i 800) 2 ] . The values of N are 
0.95 x 10~ 4 and 0.29 x 1CT 4 for the production 
angles of 18° and 5 6°, respectively. 



FIGURE 6.— Calculated differential neutron production 
cross section for 580-MeV protons incident on Be. 
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FIGURE 7.— Depth dependence of absorbed dose for 
neutrons incident on an infinite slab of tissue. 
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FIGURE 8.— Depth dependence of absorbed dose for 
neutrons incident on an infinite slab of tissue. 

The contributions of various kinds of particles 
to the absorbed dose and the dose equivalent have 
also been calculated for the case of incident neu- 
trons, and these results are given in reference 7. 

In summary , the good agreement between the 
present calculations and the measurements indicates 
that the disagreement between the previous calcu- 
lations (ref. l) and measurements was a result of 
not taking into account some of the details of the 
experiment in the calculations and not due to any 
shortcomings of the calculational method. 
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